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Abstract: 

The supersymmetry transformation relating the Konishi operator to its lowest descendant in the 10 of 
577(4) is not manifest in the Af = 1 formulation of the theory but rather uses an equation of motion. On 
the classical level one finds one operator, the unintegrated chiral superpotential. In the quantum theory 
this term receives an admixture by a second operator, the Yang- Mills part of the Lagrangian. It has long 
been debated whether this "anomalous" contribution is affected by higher loop corrections. We present 
a first principles calculation at the second non-trivial order in perturbation theory using supersymmetric 
dimensional reduction as a regulator and renormalisation by Z-factors. Singular higher loop corrections 
to the renormalisation factor of the Yang-Mills term are required if the conformal properties of two-point 
functions are to be met. These singularities take the form determined in preceding work on rather general 
grounds. Moreover, we also find non- vanishing finite terms. 

The core part of the problem is the evaluation of a four-loop two-point correlator which is accomplished 
by the Laporta algorithm. Apart from several examples of the Ti topology with two lines of non-integer 
dimension we need the first few orders in the e expansion of three master integrals. The approach is 
self-contained in that all the necessary information can be derived from the power counting finitcness of 
some integrals. 



1 Introduction 



In the M = 4 super Yang-Mills theory (JV = 4 SYM), half BPS operators are known to be protected. 
Their two- and three-point functions receive no quantum corrections other than contact terms [1, 2]. The 
simplest such operator is 

O = Tr(</>V) (1) 

where the notation refers to the action of the theory (91) in terms of Af — 1 superacids [3]. The physical 
fields are the 6 = 6 = components of the three complex chiral fields and of 

A*= \ a = D & D" (e 9V D a e- 3V ) (2) 
Kx = Va^" 1 - e^ v (D a (e^ v ^- 1 e^ v ))e-^ v , {2,3,4} 

and their complex conjugates, while the gauge field is in the [D a ,D a ] component of V. The little 
i G {1 . . . 4} is an index in the fundamental representation of the non-manifest but unbroken SU(4) 
R-symmetry of the theory. Using the equation of motion 

-\ D a D a <f> T = | e IJK [4>j, $k], I G {1, 2, 3} (3) 

of the chiral field we find 

X - D a D a O = F-AgB, (4) 
F = Tr(X^Xl) , B = Tr((e-^^e^) [0 2 , 3 ]) . 
The operator F is an SU (4) component of 

F w,ij = Tr(A ia A^) , (5) 
and similarly for B. The associated highest weight components are 

F = Tr(A«A Q ), B = Tr fr 1 [<? , <?]) . (6) 

We verify by explicit graph calculations at 0(g 2 ) in Section 3 for the highest weight components, and at 
0(g A ) in Section 5 for the 22 components F, B, that the combination 

O w = F w -AgB w (7) 

is a protected operator. Further, the contact contributions in its two-point function are related to those 
of the O case by the classical e.o.m. (3), see Section 5. 

Second, from the complex conjugate of the classical relation (3) we find 

-^D & D & K - -3gB, K = Tr (e gV ^ e~ gV (j) 1 ) . (8) 

The operator K, (the kinetic term for matter) is usually called the "Konishi operator" . 

The conformal symmetry of N — 4 SYM implies that two-point functions of properly renormalised 
operators must take the form of a power law (N is the rank of the gauge group SU(N)) 

c (a 2 N) 

(O r (x2)O r ( Xl )) = ry A( J N) , x 12 =x 1 -x 2 ; {d r (x 2 )O s {x 1 )) = : A r ^ A s (9) 
{x\ 2 ) ryg ' 

where r, s labels various operators. Quantum corrections can thus affect the normalisation c or lead to 
logarithms that sum into the "anomalous dimension" , i.e. the g 2 dependent part of A. Both equations 
only hold up to contact terms; in dimensional regularisation or related schemes this means terms of order 
e before the regulator is sent to zero. Both O and O 10 trivially satisfy the first equation in (9) because 
quantum corrections are absent, in particular A = 2, 3, respectively. On the other hand, JC must be 
renormalised. By various arguments (first from the OPE [4]) it is known that its dimension behaves like 

A^ = 2 + 3a-3a 2 + ..., a = iJL . (10) 
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K := — D h D & fC = B + ^_F+... . (11) 



Curiously, B as the "descendant" of tC under the classical equation of motion is not orthogonal to O = 
F — Ag B. Moreover, from the explicit results below we would calculate an anomalous dimension Ab = 
3 + 3 a — 9/2a 2 + . . . and hence not the same as for IC. This is in contradiction to supersymmetry: The 
N = 1 supersymmetric version of the power law in (9) is obtained by replacing x\2 by a supersymmetric 
line element . Applying the differential operator from the equation of motion on either end (thus D 2 | i D 2 \ 2 ) 
simply produces a box operator (we illustrate this in Section 5 on the simpler case O) . The dimension of 
the correct descendant is thus higher by one unit whereas its anomalous part must agree. 

According to [5, 6] in the case of K, the equation of motion is modified in the quantum theory: The 
correct descendant is 

— DaD" IC = B + -^- n 
12 g 32 7T 2 

This is an effect of the renormalisation of composite operators and not a supersymmetry anomaly. When 
JC is regularised by point splitting it is in fact possible to derive the lowest order F admixture from the 
supersymmetry variation of a Wilson line between the separated chiral fields in JC [6] . The mixture K is 
orthogonal to O at order g , and we show in the present article that it has the same anomalous dimension 
as JC up to 0(g 4 ). 

In theories without matter self-interaction the classical equation of motion would send IC to zero, 
which implies the conservation of an axial current in the Af — 1 Grassmann expansion of the operator. 
In the literature, the F admixture therefore has been coined the "Konishi anomaly" in analogy to the 
standard axial anomaly. In several models without superpotential the F term does in fact not receive 
higher loop corrections [7]. The more general case with matter self- interaction is included in a similar 
statement in [8]. 

In [9] the question was taken up again in J\f = A SYM by an analysis of the conformal properties of 
two-point functions as stated in (9). Unfortunately, point splitting is awkward to use at higher orders in 
perturbation theory. Therefore the discussion was build around supersymmetric dimensional reduction 
(SSDR) as a regulator [10], on the expense of losing the direct derivation of the "anomaly". On rather 
general grounds (the singularity structure and fractional dimension of the correlators, supplemented by 
parametric differentiation as in renormalisation group reasoning) it was shown that the "anomaly" must 
be affected by renormalisation at least in the given scheme, in apparent disagreement with [8]. 

To lend further support to the claim, we present a perturbative treatment in SSDR at the second non- 
trivial order. 1 As in [9] we resort to the conformal properties (9) to resolve the mixing. The programme 
is carried out from first principles: It is shown that the two-point function of O = F — A g B is protected 
at 0(g 2 ) and 0(g 4 ). We then impose orthogonality of 

K = Z B B + ^Lz F F (12) 

to O and the conformal form of the two-point function of K at the first two non-trivial orders in g 2 . Of 
necessity, Zp contains singular higher loop corrections. 

Upon rewriting F = O + 4 g B and rescaling, the "Konishi descendant" takes the form 

K = ZbB + ^ZoO (13) 

for which the Z-factors have been worked out in [9] in closed form in terms of the anomalous part of the 
dimension: 

Z B = cxp ^j- J dr^j , 7 (r) = 71 t + 72 t 2 + . . . (14) 



Zo = I - Z B (a, e) I W(T ) d T (j^) dr 



Lu(t) = T + UJ\T + LO2T +••• 



Note that the two Z-factors are not proportional. Apart from the anomalous dimension 7(a), graph 
calculations can only determine the finite part ui(a) of Zo- In terms of this data the poles in both 
Z-factors are given by the formulae above. In particular, Zo is singular even if = 0. 

The value oji = —3/4 was given in [9], below we establish UJ2 — 7/8. Moreover, up to 0(g 4 ) the 
singular terms in the renormalisation factors obtained by direct calculation coincide with the values 
expected from (14). 



1 Potential ambiguities in SSDR at very high orders [11] do not play a role in this work. 
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2 Loop corrections to the Konishi anomaly in SSDR 



In the next three sections of the article we work out the correlators (BB), (FF) and (BF) at the first 
three orders in perturbation theory. We use W = 1 superfields in supersymmetric Fermi-Feynman gauge. 
The regulator is supersymmetric dimensional reduction. The hard part of the problem is (BF) g s which 
contains 58 four-loop two-point supcrdiagrams. Breaking down the Grassmann algebra does not lead to 
simple numerators. We have evaluated the integrals on the computer using the Laporta algorithm, i.e. 
integration by parts. A sketch of the approach is given in Section 6. The numerator reduction has been 
independently verified by A. Pak. 

In this section we discuss the operator mixing problem between B and F drawing upon the results, 
without touching upon the details of the graph calculation. The renormalisation procedure laid out below 
is insensitive to division by regular functions of e, the decrement of the space time dimension D = 4 — 2e. 
In equation (17) we have scaled down the correlators by (BB) g o, so 



BBo 



(B(x 2 )B(x 1 )) g 



(B(x 2 )B( Xl )) g o 

etc. and we use the abbreviations (at the given orders there are no corrections subleading in N) 

g 2 N 



X = x\ 2 Ax . 



a = 



Att 2 



(15) 



(16) 



The vanishing of (BF) g i is a trivial consequence of the Feynman rules; there is no diagram. 



BB Q = 

BB 2 = 

BB 4 = 

FF = 

FF 2 = 

FF A = 

BFi = 



1 

aX' 
e 

16 g 2 ^ 

16 g 2 
16 g 2 ^- 




(17) 



[-3-3e-9C(3)e 2 + ..." 



45 



45 63 



[ 2-4e + 2e 2 +../ 



-l-12C(3)e- ( — -54C(3) ) e< 



-2 + r- 



12 + 9C(3) + yC(5) j < 2 +- 



BF 3 
BF 5 



g 



g 



[-9-3e + 24e 2 + ...] 
y + 35e+(-3 + 90C(3)) e 2 + . . . 



It follows immediately that the two-point function of O = F — 4g B is of order e at g 2 and at g A . In the 
limit e — > these are contact terms [2]. When x\ 2 ^ the one- and two- loop corrections tend to zero — 
in other words, this linear combination is protected. On the other hand, protectedness at the first two 
orders strongly point towards all-loops protectedness. Turning the argument around we might impose 
protectedness and view the absence of the pole- and finite parts of the two-point function of F — 4 g B as 
a constraint relating the e expansion of, say, (FF) to that of the other correlators. In this way the leading 
term in FF 2 and the two leading orders in FF4 must take the values given in the table. Moreover, FF 2 
must not have a simple pole and FF4 must not have a double pole. A graph calculation merely yields a 
consistency check. 

The mixture O has vanishing and hence well-defined anomalous dimension. Clearly we can construct 
a second operator K starting on B that will have to be orthogonal to O if it has non- vanishing anomalous 
dimension. To lowest order in g 

gN „ 



K = B + 



32 7T 2 



(18) 
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satisfies this constraint, but we also see that 

(K(2)0(l)) gl = (B(2)B(l)) g o4gelog(x 2 12 fl 2 x ] 



(19) 



from which it is clear that we cannot expect the orthogonality constraint to hold beyond O(e ). The bare 
two-point function of K is divergent at 0(g 2 ) because the pole in BB 2 cannot be compensated by the 
finite contribution FF . Wc therefore renormalise as 



K 



Z B B 



gN 

32 7T 2 



Z F F. 



(20) 



Z B = 1 + a 



bn 



+ a' 



&22 &21 



1 + a 



+ /lO + OL 



+ 



/22 



hi 

e 



20 



Finite terms in Zb could be absorbed into an overall rescaling and a corresponding change in Zp\ omitting 
these amounts to fixing the normalisation. 



(K(2)K(l)) g2 
(B(2)B(l)) g0 



2a — + BB 2 + 



( gN 



V32 7T 2 



FF 



(21) 



= a 



bn 



3 -31og(X) 



0(e) 



On the other hand, in the renormalised QFT conformal invariance implies that the two-point function 
has the functional form 



(K(2)K(l)) g2 
(K(2)K(l)) g0 



I + aai + a 2 a 2 + . . . 



1 + a( ai - 7i log(X)) + 0(a 2 ). 



(22) 



In the first equation we can put e to zero after adjusting bn. It is then possible to equate the last two 
lines. We learn 

7i = 2&H = 3, ai = --. (23) 



Next, let us put 



n X e 

BF 3 = g [B 31 -3e + 0(e 2 )] 



(24) 



for the moment and consider the orthogonality constraint at 0(g 3 ). Using BB = 1, BFi = we obtain 



(K(2)Q(l)) g3 
(B(2)B(l)) g o 



= BF 3 -4gBB 2 -4ga— + ^J^FF 2 r 



32 it 2 



gN 

32 7T 2 



■ a 



fi 



+ ho FF 



(25) 



= ga 



- (B 31 + 6 + 4/ n ) + log(X) (B 31 + 12 - 4/u) + (7 - 8 f n + 4/ 10 ) + 0(e) 



This vanishes up to O(e) if 



B31 = -9, 



11 



10 



(26) 



The singular term /n/e in Z^? is necessary, because we need to solve two different equations in order 
to eliminate the pole and the logarithm. This fixes both the fn term and the pole term in BF 3 . 
According to (17) the graph calculation does indeed meet the requirement B 3 \ — —9. The finite pieces 
in BF 3 , BB 2 , FF 2 and fn times the subleading order in FFo determine the remaining coefficient /10 = 
— 1/4 7^ 0. Hence in this scheme the Konishi anomaly is affected by renormalisation. 
Up to this point we reviewed the analysis of [9]. Next, we look at 



(K(2)K(l)) g , 
(B(2)B(l)) g0 



2bn 1 
BB 4 + a — - BB 2 + a 2 



b 2 11 +2b 22 2 b 2 i 



+ 2^BF 3 
32 7r 3 



gN 

32 7T 2 



FF 2 + 2 a 



e 

hi 



(27) 



+ ho FF ( 
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or by substituting the explicit formulae 



(K(2)K(l)) g , 
(B(2)B(l)) g0 



where we have put 



a 



1 



27 



1 



-= \B i2 + 2b 22 - -r + - log(X) (2B 42 - 9) 



BBa 



2&21 + 



aX ( 



+ log(X) 



45 



21 



17 
Y 



■C(3) 



45 63 



B i2 + -e + - + -C(3) 



(log(^)) 2 

+ 0(6) 



+ ... 



2B 42 - 2 I 4 



(28) 



(29) 



for now. The elimination of the singular terms leads to three conditions, one of which we solve for B42. 
Upon equating with (22) we find 



B. 



42 



^22 



72 



»21 



«2 



17 9 



(30) 



Once again, the leading coefficient in BB 4 stated in (17) does in fact take the right value. 

It remains to analyse the orthogonality constraint at 0(g 5 ). The system is overdetermined as before 
so that we start with 



BF* 



aX' 



[ J B 5 2 + 5 51 e+(-3 + 90C(3)) e 2 + •••] . 



(31) 



(The pole part of FF4 is already constrained by the protectedness of O 10 ). We put in the other correlators 
from (17) and the coefficients of the Z-factors already derived. 



(K{2)Q(l)) g , 
{B(2)B(l)) g o 



ga 



-2 (B52 + 4/22 - 18) + - log(X) [2B 52 - 4/22 - — + 



03 



(32) 



135 



1 



+ (log(X)) 2 {2 B 52 + 2 f 22 - —j + - {B51 +4 /21 - 8 / 22 - 30) + 
+ log(X) (2 B B1 - 4 /21 + 8 / 22 - 75) + (4 f 20 - 8 f 21 + 4 / 22 - 7) + 0(e) 



The correlator vanishes up to 0(e) if 2 
33 

-B52 = 7 B 51 = 35 , 



h 



f: 



20 



(33) 



In conclusion, all three 0(a 2 ) mixing coefficients in Zp are no n- vanishing. Interestingly, £(3) cancels 
from f 2 Q. The fact that the sum of the graphs in (BF) g s reproduces B B2 , B$i is extremely non-trivial by 
looking at the orders 1/e 3 and 1/e 2 in the explicit results in the equations (55) and (60). The constraints 
imposed by conformal invariance after renormalisation give an excellent test of the graph calculation. 

To make touch with the discussion in [9], we write F = O + 4g B and divide K by l+a/2 — a 2 / 8 
in order to eliminate finite contributions from the shifted renormalisation factor Zb + aZp/2 which now 
multiplies B. The renormalised operator mixture fits the general form 



K 

Z B 
Zo 

predicted by (14), with 



Zr B 



SN ZoO 



32 7T 



(34) 



7i 



1 + a 
l + a( 



2e 
4e 



a 



8e 2 



72 
4e 



+ 0J1 



)+a 2 



7i 
24 e 2 



7 2 + c^iTi 
6e 



w 2 



71 = 3, 



72 = -3 , 



uji = -- 



L0 2 = 



7 



(35) 



Our effort shows that the renormalisation scheme developed in [9] is operational also at the second order in 
a, and it yields the previously unknown number u 2 which could be called the two-loop Konishi anomaly. 



2 We could have determined 622 from here, too. 
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3 (BF) 



The complete set of superdiagrams is displayed in Figure 1 and Figure 6. The diagrams in Figure 1 are 
special in that they have two cubic non-abelian vertices or one non-abelian four-vertex. They all turn 
out to be derived topologies 3 of the first diagram Go so that they can conveniently be summed into one 
effective numerator Go- 

Go 




Figure 1 

Suppose that we start by Wick contracting only the matter part of the diagrams, i.e. the operator 
B(xy) on the cubic chiral vertex at point 6 and the respective matter/YM vertices at point 4 and 5, 
where present. There are three "matter parts": 



Mi 



Mo 



Ma 




V[A) 



V(5) 




6 V(5) 

Figure 2 




The matter line has an arrow pointing from the antichiral to the chiral end. One way of writing this 
propagator is 



mm) 

Here </>(«) = (p(xi,8i,§i), and 



D.D D.D\iIlij , 



n, 



D.D 



Wj) 



1 



D a D n 



D.D 



1 



DnD a 



(36) 



(37) 



-IF; 



where Xij,9ij0ij denote coordinate differences. The Yang-Mills propagator is (V(i) V(j)) 

For this class of graphs we will use "Z?-algebra" to derive the numerators. The tools are partial 
integration, the commutation relation {D a ,Da} = —2id a a and the shift rule 



Di Tly 



(38) 



The idea is to free all but one line of each loop from spinor derivatives. Due to the shift identity and 
the possibility of reordering via the commutation relations, we may reduce to spacetime derivatives and 
maximally two chiral and two antichiral spinor derivatives on each 11^ . On the last line in each loop we 
need to keep only terms with all four derivatives, since 



(D.DD.D S(6ij) 5(§ij)) = 8(e i3 )5(e i3 ) 



(39) 



3 topologies that arise by cancelling a line 
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while less than four derivatives between two (5-functions vanish. 

Let us evaluate matter part Mi with this technique: We absorb the D.D of the 7-6 line into the 
measure at the chiral vertex 6 and partially integrate the D.D away from it onto the 4-6 and the 5-6 
lines. We then shift all spinor derivatives to points 4 and 5 respectively. Last, we partially integrate the 
spinor derivatives away from the 4-6 line, onto V(4) and the 7-4 line, likewise at point 5. The result is 



Mi 



-6N 2 g 3 



Tr [(D.DV (4)) V (5)} D i6 D 57 + Tr [(D.D V(4)) (D & D a V(5))] D 46 - d« 7 a 
Tr[(D.DV(4))(D.DD.DV(5))] n 46 + 

Tr [V (4) (D.DV (5))} D 56 U„ + Tr [(D & D a V(4))(D.D V(5))] a 56 l -df 7 a 
Tr [(D.D D.DV(4))(D.DV(5))] D 56 - 
Tr[(D a V(4))(D a V(5))} ^D 46 a 5e - 

Tr [(D a V(4)) (D.D D & V(b))} U i6 % - 8^ + 

Tr [(D.D D & V(4)) (D a V(5))] D 56 % - 8™ - 



Tr 



(D.D D & V (4)) (D.D D? 1/(5))] - 



(40) 



In the last formula dij only acts on the i — j line, at point i. For notational convenience we have omitted 
n4 6 n47ll 56 n57n 6 7 and the integrations. Similarly, 



Mi 
67V V 



= Tr [V (4) (D.DV (5))} D 56 D 47 + Tr [(D & D a V(4))(D.D V(5))] D 56 % - 0* Q - 



Tr [(D.D D.DV(4))(D.DV(h))] □ 



56 



Ti[(D.DV(4))V(5)] D56D45 - Tr \(D a V(4)) (D a V(5))] ~n 56 D 45 



Tr [(D.D D A V(4)) (D a V(5))] □«, - 8[ 



a. a: 

'45 



(41) 



Next, D 56 leads to the same derived topology in Mi and M 2 (by shrinking the 5-6 line to a point) so 
that the "numerators" can simply be added: 



Mi + M 2 + (M 2 : (4 o 5)) 
6N V 



Tr 
1 



1 



Aiind(4)A£ n (5)j 2 9 46 Q ^ + (Tr[(D.DV(4))V(5)} (42) 
Tr [(D a V(4)) (D a V(5))} + Tr [V(4) (D.D V(b))]) U i6 D 56 

D.DD & V. On the 



with the linear part of the antichiral physical fcrmion in the gauge multiplet Xf in 
other hand: 

" U:i = -D.DTr[V(4)V(4)} n 4 6 



(43) 



QN 2 g 3 

Due to Greens' function equation n 56 = +5 8 (z5e) this cancels the C^Dse terms in (42). In Figure 
1 we can therefore work with one effective matter part 

M = 37VVTr[A lind (4)A^ n (5) 



a 46 °5&af} ■ 



(44) 



Note that putting 9 7 — 9 7 — also selects the 6 — 8 — component of Ai; n at both open ends. 
To show consistency with [9] we recompute the leading 0(g 3 ) contribution to (BF). With 



F = Tr [A a X a ] A° 



■D.De gV D a e- gV 



we find 



(B(7) F(l)) g3 = -24 g 3 N 2 (N 2 1) d% d 56a$ d?£ d li0& 



(45) 



Once again, for notational convenience we only gave the numerator and omitted the integral itself. 
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Throughout the paper we employ supersymmetric Fermi-Feynman gauge. The last result would in 
fact have been immediate in Wess-Zumino (WZ) Fermi-Feynman gauge. However, the D-algebra can be 
done with the aid of a computer programme so that the non-supersymmetric gauge does not offer any 
advantage for the g 5 part of the correlator. 




1 Jf-"" 






1 4 








' 5 


3 ~- 





Figure 3 

The first picture in Figure 3 displays the effective matter part Mi contracted onto F(l), or equivalently 
the one graph in (B F) g a in WZ gauge; here the numbers label the spacetime points. The second picture 
reproduces the momenta assignment from the discussion of three- loop p-integrals in [12], where this 
topology is called BU. Fourier transform leads to the identification 



d 



14 



ip 6 . 



di5 -> -i P7 , d 4 e -> -ip4 ■ 



^56 -> -ips ■ 



(46) 



The negative sign of the five matter propagators accounts for an extra minus, so that in momentum space 



(BF) g3 = 



-48 g 3 N 2 (N 2 
12 g 3 N 2 (N 2 



!) [(P4-P5)(P6-P7) 

■1) 



(P4-P7){P5-Pe) + (P4-P6)(P5-P7)] 



(47) 



(47T) 



1 28 1 166 n . u 

— + -5 1 5~~ + O e 1 

3 e 3 



o 2 



1 — ) 



The second line shows the result of Mincer in its MS 1 mode. 4 

After this digression let us return to the diagrams of Figure 1. Since combinatorics is linear, the 
collection of ten supergraphs can be viewed as the effective matter part contracted onto the pure Yang- 
Mills correlator 

Yi = (A^ lin (4)A^(5)F(l)) ff2 (48) 
where a, b are colour indices. The latter contains the following graphs: 




Y 2 

A(4) 



AM) 




F(l) 



Figure 4 



Y 3 

A(4) 



A(5) 




F(l) 



Y 4 

A(4) 'Hn 
A (5) ^ 



F(l) 



The D-algebra for this sector is cumbersome; as an example we sketch a strategy for the first diagram 
Y\ . There are three derivatives on each outer leg and four at each vertex. We resorted to a Mathematica 
script for Wick contraction and superspace algebra. One may choose to shift the derivatives to point 2 
and 3, respectively, and then clear the 1-2 and 1-3 lines by partial integration. The derivatives on the 2-3 
line are subsequently shifted to one end and reordered. Due to formula (39) only terms with D.DD.D 
need to be kept. Next, we regroup the derivatives on the 2-4 and the 3-5 legs. Upon putting the outer 
theta variables to zero only terms with four spinor derivatives on both the 2-4 and the 3-5 legs survive. 

In a similar way, 



Y 2 



2Ng 2 6 ab 

Y± 
2Ng 2 S ab 



~ $ (Di - □« - Dm) + (d?f d 2 4p* - dH 9 15/36 )) - 3 5% □ 
Si (pi + H D 1S + 1 D 24 ) + {d^ 9 25ffA - dH d 24fj& ) 



(49) 



where d\ is a total derivative w.r.t. x\. Clearly Y 2 , Y3, I4 arise from Y\ by [H13, D^, ^23, respectively. 
Here the Green's function equation brings in an additional minus sign because (V(i)V(J)) = — n^-. Upon 



The definition of the mass-scale fi is given before equation (55). 
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adding these contributions into Y\\ 



Yx 



2Ng 2 5 ab 



4 n^djf 9 3 5/3d - dH d 2i06l ) - (D 24 + □ 35 )(9f 13 d 2306l - 8^ d 106l ) 
-4Di(Di - D 24 - Das) - (D 2 4 - a 35 ) 2 + (D 12 D 35 + D^f^) 
+ (□12024 + D 13 D 35 ) - + 035)023 



(50) 



As another cross-check we derive (F(A) F(l)) g 2 from these formulae: 

2g 2 N(N 2 -l) [-8D? + 7ni(n 2 4 + n 3 4) 



<5|<5 ab Yi(5^4) 



3 (□i 2 D 3 4 + Oianad)] 



(51) 



When identifying the points tadpole terms were put to zero. Finally, the mirror image of Y 2 , Y 3 arising 
from the cubic part of F(4)| g and the tree-like contribution from (F(4)| s F(l)| g ) g 2 have to be added. The 
result is 

<F(4) F(l)) g2 = 32(N 2 - 1) g 2 N (D 12 D 3i - \ U\) (52) 

consistent with the results in [9] for the SU(4) component F = Tr((V"$ 1 )(V a 4> 1 )) of F 10 . According to 
(63) in the next section, the second term in the last formula is a contact contribution proportional to £(3). 
The two box operators in the first term break the integrations and we find —32g 2 N(N 2 — 1) (0(4)</>(l)) 3 
due to the opposite sign of the matter propagator. In the combination F — AgB this is cancelled by 
16 g 2 (B(A) B(l)) g o as required by protectedness. 

By putting formulae (44) and (50) together we obtain the desired effective numerator G for the first 
topology in Figure 1. Graph G is of topology C in the list in Figure 5 (the numbers in the figure label 
the momenta p\ . . -Pio)- 



B 



D 







Figure 5 

In Section 6 wc comment on the realisation of the Laporta algorithm [13] with which we evaluated 
such four-loop diagrams with rather general six-derivative numerators. In topology C we have chosen 
the momenta q,Pi,p 2 ,p 3 ,P4 as independent. To calculate Gq we identify 



di 



iq, 



d 2i ipi , 



0; 



2.1 



<9 46 -+ i P3 ■ 



a 



5(i 



IPA ■ 



(53) 



Between these one may form ten dot products. In total there are ten interior lines, so that six of the dot 
products may be eliminated in favour of the four others and the squares of the dependent momenta. In 
terms of p 2 .p 3 , p 2 -P4, P3-q, Pi-q and p\ . . .p 2 Q, q 2 the effective numerator has 47 terms. Flip symmetries 
of the original topology and derived diagrams could be used to simplify. 



G = 12g 5 N 3 {N 2 -l) [8q 2 p 3 .qp 4 .q - 2 (p 2 .p 4 p 3 .q - p 2 .p 3 p 4 .q)(p 2 + p 2 ) + (42 terms) ] 



(54) 



(The overall negative sign from the odd number of matter propagators has been taken into account.) The 
omitted terms all have at least two distinct complete squares. Each pf, i = 1 ... 10 cancels a line in the 
denominator and thus leads to a simplification within the Laporta approach (although pf is generally not 
better than p\). The first term in the square brackets is thus the hardest to analyse; an algorithm that 
solves this problem can also calculate the other pieces. 

In the MS scheme fractional powers of 4 and it and the Euler constant are absorbed into the mass 
scale by putting p 2 = fi 2 Aire~ lE . One may further scale away £(2) by the redefinition fi 2e = /2 2e (l — 
e 2 7r 2 /12 + 0(e 4 )). In this convention our result for the sum of the diagrams in Figure 1 takes the form 



Go 



12. g 5 ^ 3 ^ 2 - 1) 
(4tt) 8 



19 j_ 

72^3 



53 J_ 

16^2 



/7957 53 



V 288 



C(3) 



O(e 



-4e 



(55) 



At the first glimpse not all diagrams in Figure 1 can in fact be drawn into topology C. We have shown 
above that the numerators of the matter parts M 2 , M 3 have d'Alembert operators acting onto the 5-6 
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and 4-6 lines, respectively, common to all terms. These are marked by grey boxes in the graphs of Figure 
1. Contracting the corresponding lines yields topologies that can also be derived from C. 

We have indicated such "topology changing" box operators on the diagrams of Figure 6, too. For 
instance, graphs G31 and G32 (the one- loop correction to the chiral vertex in matter part Mi) have 
originally neither of the four topologies A, B, C, D. The grey box operators are positioned such that 
G31 can be obtained from B by a factor p\, and G32 from A by p 2 or equivalcntly from D by p\. 

Figure 6 contains several classes of diagrams that could be summed prior to evaluation like the graphs 
of Figure 1. The graphs Gi 6 . . . G 20 , G 2 i . . . G 25 , G 26 • • • G 30 and G 3i , G 32 ; G 33 , G 34 and finally G 40 , G41 
are all vertex corrections relating to the matter parts Mi,M 2 ,M 3 contracted onto F(l). Even in the 
most complicated graphs Gi6, G19, G 2 o (type D) there is at least one square of an interior momentum in 
each numerator term, hence the complexity is lower than for Go- There is no calculational disadvantage 
in individually treating the diagrams, so that we separately list the contribution of every supergraph in 
order to make the whole calculation more verifiable. 

The longest numerators for the A, B topologies arise from the graphs G3 and G10 in Figure 6 with 
134 and 183 numerator terms in our momenta assignments, respectively. These as well as the numerators 
of Gi, Gq, G-j contain terms with no squared momentum. Yet, in all such terms at least two of the 
three non-trivial dot products involve the outer momentum q. The numerator reduction can then still be 
accomplished with a slightly restricted matrix; terms with one or no q are not needed in the ansatz since 
q cannot be made to disappear by differentiation w.r.t. loop momenta. 

In our way of organising the Laporta algorithm by stepwise elimination of lines (see Section 6), 
topology A is the hardest case because it contains only the triangle between p4,p 8 ,p 9 . Using the triangle 
rule (84) we immediately fall upon two derived topologies that must be further reduced by the Laporta 
ansatz. In topology B we use (84) on the triangle consisting of p4,p-j,pg,, which produces a new triangle 
in every term. The step to the eight-propagator level is therefore more direct. Topologies C and D both 
contain two disjoint triangles. 

Graph G4 is in the same way a subgraph of G3 as the matter part M 2 belongs to Mi; likewise, G19 
and G 2 g go together, and G10 and G 24 . On the other hand, G 24 naturally belongs to a class of vertex 
corrections, while G10 does not. It is thus not clear how to assemble all the superdiagrams in Figure 6 
into well-defined classes; likewise Gn and its subcase Gi 2 could be attributed to cither of the three pairs 
with one non-abelian vertex. 

The "D-algebra" technique is not convenient in diagrams with too many matter lines, because ac- 
cording to (36) every such propagator brings in four spinor derivatives. In the diagrams of Figure 6 we 
have rather evaluated the Yang-Mills part (in Gi and similar diagrams the two lines attached to F(l), in 
G3 and G11 etc. the complete gluon part) by D-algebra and then taken the remaining spinor derivatives 
to obtain the Grassmann expansion at 6\ = 61 = 0. The alternative form of the matter propagator 



(4>{i) <f)(j)) = e - l (( e ^^+(M^)-2(M*,e»)) 



1 



-47r 2 4 



(56) 



easily lends itself to Grassmann expansion. The complete numerators are found by taking the product 
with the expansion of the Yang-Mills parts and integrating out the spinor variables. 5 
Finally, from the sum of Go and Gi . . . G 4 s of Figure 6 / equation (60) we find 



(B(7)F(l)) g5 = 



12 g 5 N d (N 2 - 1) 



11 1 899 1 / 15043 nn „ n .\ 1 nl ' 



in momentum space. The backward Fourier transform 6 

1 1 r(f-a) 



x e 



(2n)D J (g2)a ^ T{a) ^§-a 

yields (// 2e = (^ 2 7re 7E ) e (l + e 2 7r 2 /12 + 0(e 4 )) in configuration space) 



<S(7)F(l)) fl . 



g 5 7V 3 (jV 2 - 1) 

(4^2)5 



33 70 

^ + — + (-6 + 180C(3)) + O(e 1 ) 



.5e 



(57) 



(58) 



(59) 



5 We have computed Go also by diagram-wise evaluation using this technique. 

6 The forward transform lacks the factor 1/(271-)-° and has the roles of q,x and the sign in the exponent exchanged. 
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1 ( - r 2 <J"3 ^4 




Figure 6 

Diagrams with vanishing colour factors (w.r.t. SU(N) gauge group) or one-loop propagator bubbles 
were omitted. Graphs G22, G23, G27, G28, G33 . . . G48 do not need numerator reduction by the Laporta 
ansatz, though the algorithm is needed here, too, to determine some T 1 configurations (see below) with 
one or two non-integer exponents. 

In (60) we give the momentum space result for each graph including all signs and factors. The e 
expansions in the table must be multiplied by a global factor — 12 g 5 N 3 (N 2 — l)/(47r) 8 * q 2 [q 2 /jl 2 ^ e . 
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4 (BB) g , 

At tree-level 



2N{N 2 - 1) 1 



W (s? 2 ) 



_> , 3 (^/';) :; ' (l + 0(e 3 )) 



(61) 



where the e 3 C(3) + • • • correction arises from T(l — e) in the numerator of the modified propagator in 
dimensional regularisation. At order g 2 there are two graphs: 



(S(2)S(l)) ff2 = 12 g N (N — 1) (</>(2)0(l))„o (ff a + H 2 ) 



(62) 



Hi 



Figure 7 

Hi 

whereby in configuration space 
(B(2)B(l)) g2 = 




1 



(4tt) 4 
1 



14 

+ - + 12 + e 



(63) 



, 4 [3C(3) + 0(e)] . .., 



,g 2 iV 
"4^2"' 



B(2)2?(l)) 9 o 



- + 3 + 9C(3)e + 0(e 2 ) 



(64) 



The diagram B% is of topology Ti in the nomenclature of [12] . Its numerator is just an outer box operator 
so that we can take over the result (87) elaborated below. The leading term is of order e° in momentum 
space, which means 0(e) in configuration space. It is therefore a contact term [2]. Note that the contact 
contribution in (52) is exactly H% with a second outer box operator, which does not touch upon the order 
in e of the leading term. 

On the other hand, w.r.t. the 0(g 4 ) part of {BE) we are only interested in the singular part in 
momentum space, or equivalently the singular and finite pieces in configuration space. On comparing 
to a certain protected correlator we can spare most of the work: Expanding in N , protected correlation 
functions contain several linear combinations of graphs in which the x space singular and finite terms 
must cancel. One such sum of graphs is also present in {BB) g i, with identical relative coefficients for 
the relevant graphs. Using the work of [14]: 



<<K2)<Mi)} 9 o 



12 g 4 N 3 (N 2 - l r 
where Si contains the six four-loop graphs H3 



B(2)B{l)) g , = (0(2)0(l)) 5 oS 1 



H 



(65) 



■#8 




and 



£3 = 



Figure 8 

(P(2)P(1)V 



P = Tr(0VW) + -Tr(0VW) 



(66) 



is a sum of graphs in which only a contact contribution survives because P is an SU (4) component of 
the half BPS operator Tr^ 1 (j) 1 (f) 1 (j) 1 ) . Hence in x space 



{B(2)B{l)) g 



12 g 4 N d {N 2 -l)Ei + 0(e) . 



(67) 
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The first five graphs in Si are simple to evaluate with the methods developed for (BF) g s . For convenience 
we used Mincer for the three-loop ladder subgraph in H^. Including all signs and factors 
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to be multiplied by an overall 12 g i N z (N 2 - 1)/(4tt) 8 * <? 2 (g 2 //i 2 
space 



Summing up and translating to x 



(B(2)B(l)) g * = 



4tt 2 



(S(2)5(l))„o 



9 1 45 1 /45 63 .,„, 

2^ + ~7 + T + T c(3) 



O(e) 



(*i 2 ^) 2£ 



(69) 



5 (CC) g4 and {FF) g 4 

Let us first consider the two-point function of the half BPS operator 

O = Tr(</>V) . 

The two leading orders are 

((5(2)0(1)) = 2(iV 2 -l) [W2)^(l)) 2 o + 4g 2 NH 2 + 0(g i )] 
while the 0(<? 4 ) part receives contributions by the following graphs: 
h h h h h 



(70) 









Figure 9 

The blob in J5 means the one-loop correction to the cubic vertex g J Tr^^V, $/]) and the blob in Iq 



denotes the two-loop correction to the matter propagator. 



h 
h 
h 
h 
h 
h 



C(3)' 
C(3)' 



"4C(3)i 
2C(3)i 



h + 6 f( 3 ) 



HO 



15 



60 

+ 6C(3) + 
20C(5) 
-30C(5) 
- 24C(3) + 

M + 12C(3) 



5C(8) 



30C(5) 



■0(e) 
0(e) 

■0(e) 

■0(e) 
0(e) 

■0(e) 



(71) 



The e expansions in the table come with an overall factor 8g 4 N 2 (N 2 — l)/(47r) 6 * (q 2 / fi 2 )~ 3e . The sum 
of the lines in (71) is 25£(5) + O(e); hence only a contact term remains as required by protectedness. 

We only give the sum of the propagator and vertex corrections because the loop order here is not 
higher than three so that we could use the Mincer programme to compute the graphs. Moreover, the 
vanishing of the simple pole in the sum has been checked in [2]. 

Note that I\ . . . I3 are only of order e _1 ; in a; space they are therefore individually finite, while J4 is by 
itself a contact term. Likewise, the five graphs contributing to J5 are all separately finite in x space. In 
N = 4 SYM non-protected operators have to be renormalised only due to diagrams with subdivergences 
coming from the outer ends; as typical examples we mention the pictures of H\ or H3, H$ above. For 
protected operators these divergences are absent. Due to the improved power counting in superspace, 
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connected graphs tend to be finite when they have sufficiently many outer legs. On the other hand, the 
individual graphs contributing to the (overall finite) two-loop matter blob have poles up to second order 
in momentum space which the "kinematical" third loop in Iq even promotes to 1/e 3 . 
Collecting terms, in configuration space 



1 a (6C(3) e + 0(£ 2 )) {xl^lf 

+ a 2 ( 7 4a5)e + 0(^)) (xl 2 tir + 0(a 3 ) 



(0(2)0(1)) = (0(2)0(l)) 9 c 



Next, writing = Tr [(e gV (f^e^) (e sV '</> 1 e 9V ')] and using the equation of motion 



(72) 



g 2 N 
4^ 



D.D (e 



) = 9 [<^2 



(73) 



we derive 



- D a D a O 

F 
B 



= F - AgB, 



= Tr [(D a (e-^^e^)) (D a (e"» v cfe^))] , 



Tr[(e-^^e^) 



Recall that the F involving two Yang-Mills fcrmions which we used before and the representative F 
arc different 5J7(4) components of the same operator, likewise for B and its TV" = 1 chiral companion 
B = Tr(<W,0 3 ])- 

The Grassmann expansion of the correlator (0(2) 0(1)) is given by the same exponential shift operator 
as for the superspace matter propagator: 



(0(2) 0(1)) = e K( () ^ 1 S 1 )+(e2d 1 e 2 )-2(e 1 d 1 § 2 )) {l)) ei 2=§1 2=0 

It follows that 

((F - AgB)(l)(P - Ag$)(2)) = AD.D\i D.D\ 2 (0(2) 0(1)) = -4 Di (0(2) 0(1)) . 
In particular, for the 8 = 8 = component we should find to leading order in e 
(OF - 4gB)(l)(P - AgB)(2)) 



(74) 
(75) 
(76) 



64(^ 2 - 1) (x 2 - u 2y 
(4tt 2 ) 2 (x 2 12 ) 



75 

(1 + 0(e)) -ae6 C(3) (x 2 2 A 2 ) £ + a 2 ^ C(5) (x 2 2 /i 2 ) 



The tree-level part is obviously right for both F and F, whereas the O(a) contact term agrees with (52) 
because F and F have interchangeable two-point functions by R-symmctry invariance. We now wish to 
verify the 0(a 2 ) contribution in the last equation by a direct graph calculation. We prefer F for this 
purpose because F = X a X a + . . . would, of course, heavily involve the Yang-Mills sector, while we have 
seen above that the D-algebra creates quite some work already at 0(a). 
Upon expanding the exponentials in the definition in (4) 



Tr(Vy A^ 1 -2g[D a V, ^]D a 4> 1 + 

g 2 [D"V, 1 ] [D a V, 1 ] + g 2 [[D a V, V},^]D a ^ + 0(g 3 )) . 



(77) 



In the last formula D a acts only on the field immediately after it. In other words, in the given frame the 
chiral spinor derivative becomes Yang-Mills covariantized: 



V" = D a - g[D a V, .] + y -[[D a V, V], •] + 0(g 3 ) 



(78) 



(On the other hand, the constraint D (f) 1 = is not modified.) 

The top graphs in (FF) gi g = § =0 are like in Figure 9 but with a partial spinor derivative on each 
outer leg. These graphs have the e expansions I\ . . . Iq in the table below, which are visibly not overall 
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derivatives of h . . . Je . Some new subdivergences have to be compensated by diagrams arising from the 
higher terms in (77). Putting 9 = 9 = at the outer points removes supergraphs with more than two 
outer V fields. Further, most ways of placing D, D on the fields at the outer points lead to vanishing 
results. The remaining extra diagrams are given in Figure 10. We have put the derivatives onto the 
graphs where their position could be ambiguous. 



h h h ha hi 
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Figure 10 
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In the table we have omitted the overall factor F = 64g 4 N 2 (N 2 — l)/(47r) 6 * q 2 (q 2 /p. 2 ) e ■ The sum 
of h ■ ■ ■ I19 yields the first line in the next equation. 



(FF) gl 


= F 


1 1 


4 

19 4 

e 


-3C(3) 


or. 

+ yC(5)+0( £ ) 


-8g(FB) g s 


= F 


3 1 

+ 27 2 


+ H + 83 
e 


+ 0(e) 






16 g 2 (B B) g 2 


= F Q 


1 

e 2 


10 -64- 
e 


3 cm- 







17 



Upon adding up 

((F - AgB)(F - AgB)) gi = F ^((5) + 0(e) (81) 

which reproduces the 0(a 2 ) term in (76) after Fourier transform. 

In conclusion, at the given orders the operator F — 4g B does not need renormalisation. Its anomalous 
dimension is zero at one and two loops. At loop level the two-point function does pick up contact 
contributions; remarkably the O(o) contact term has normalisation proportional to £(3) while the 0(a 2 ) 
term comes with £(5). The equation of motion relating the operator O = Tr^ 1 ^ 1 ) to F — 4gB is 
apparently purely classical: To leading order in e the normalisation of the contact part is compatible 
with superspace differentiation. 



6 Integration by parts (IBP) in dimensional regularisation 



We follow the last reference in [12] in exposing the fundamental idea. Suppose that any given diagram 
has a triangle subgraph like the first picture in Figure 11. 
p p p 

1 2 




Figure 11 



The integral associated to the subdiagram is 

I(ao,Pi,p 2 ,ai,a 2 ) = [ d D p - a - a . (82) 

J (p 2 ) ao [ip+Piff 1 (A) ((p+P2) 2 f 2 (A) 

Discarding boundary terms 



/ 



d D p d p 5 — s = (83) 

>T° {{P+Pi?t (Pi) ai {(j>+P2?t iA) a2 



because the integrand is a total derivative. By working out the differentiation: 

I(a ,p 1 ,p 2 ,a 1 ,a 2 )(D-2a -(3 1 -(3 2 ) = (84) 
Pi (I(a - l,/3i + l,fi 2 ,a\,a 2 ) - I(a ,(3 1 + l,/3 2 ,ai - l,^)) + 
I3 2 (I(a - l,/?i, A> + l,ai,a 2 ) - J(a , /?2 + l,a 1: a 2 - 1)) 

The generalisation of the formula for polynomial numerators p" 1 . . .p Vn is straightforward. In every term 
on the r.h.s. of (84) the exponent of one of the propagators of the triangle's top line is diminished by 
one. To illustrate the use of the formula we compute the integral Ti (without numerator) given in the 
second picture of Figure 11. Using (84) on the left triangle in the diagram produces an overall factor 
1/e times the third minus the fourth picture. Our convention is that a line without any extra symbol 
has exponent equal to one. The triangle rule can remove such a line, but in exchange it augments the 
weight of one of the lines by one. This is customarily denoted by a dot. The resulting new diagrams are 
both elementary: One is a convolution, the other the product of two one-loop integrals. The elementary 
one-loop building block is 

where 

G(a,p) = (4nY R(a) R(p) R(D - a - p) , R(a) = T< ^ D / 2 ~ a "> , ( 86 ) 

T(a) 

For the generalisation to integrals with numerator polynomials we refer to [12]. It follows 



\ G(l, 1) (G(2, 1 + e) - G(2, 1)) = ^ ( 6 C (3) + + 12 C (3)) e + . . .) 1 ^ 



-it 
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in the MS convention explained above. Almost all three-loop "p-integrals" (propagator type) can be 
calculated by this trick [12]: By way of example, to solve the BU topology displayed in Figure 3 one may 
start with one of the visible triangles and then iterate the procedure on the resulting Ti (sub)diagrams. 

Let us now turn to the four-loop topology A of our (BF) g & problem. We can use the rule of the 
triangle once: 




Figure 12 



The two new topologies contain no further triangle. Nonetheless, the IBP technique can be pushed 
further by the Laporta ansatz [13]: The integral in the middle of Figure 12 now has propagators with 
power one for, say, p\ to p s and of power two for pg. The most complicated numerators we encounter 
have three non-trivial dot products. We write the complete set of identities 

/n n P^ num 

d D Pl . . . d D p 4 d p » 3 = (88) 

pi ■■■pi (pi) 

where i G {1 . . .4} and pj can be one of the loop momenta or the outer momentum q. Between the 
four loop momenta and q one can form fifteen independent square invariants. Let the first ten be the 
squares of all the interior momenta and of q, then we have to choose five further mixed dot products. 
Mixed products containing q bring an advantage. We write the 20 identities (88) for any such numerator 
num built from three square invariants that does not by itself cancel a line in the denominator (thus 
it can contain the mixed dot products, q 2 or one power of p|). Next the differentiation is worked out 
just as in the case of the triangle rule. The result is a large homogeneous linear system for a basis of 
integrals, which one may reduce by Gaussian eliminination. The matrix is initially very sparse and the 
elimination has the surprising property that the relative order in e between the terms in the same line 
remains relatively stable throughout the steps of the algorithm. 

Due to memory limitations (1 GB on a power PC and later a Xeon, of which we needed only about 
one half) we organised the task in a recursive way: In the first step, except for the integral with a sought 
numerator, all nine-propagator structures are eliminated. We obtain a linear equation relating it to cases 
with eight or less propagators (the cancellation of a line works as in the case of the triangle rule). The 
set of eight-propagator configurations is given in Figure 13. The differentiation in (88) produces a second 
dot, which we did not indicate on the pictures because it can be placed on any of the interior lines. 
By dimensionality the eight-propagator integrals still have numerators with maximally three kinimatic 
invariants, or less if there is a further cancellation with a denominator term. 

M 35 M 36 




Figure 13 



On this level, M 35 and M 36 once again need numerator reduction by the Laporta ansatz. There is a 
new feature arising here: Apart from the desired numerator term one can eliminate all eight-propagator 
cases but one. The conventional choice for the remaining integral is M 35 or M 36 , respectively, with trivial 
numerator q 2 and all exponents equal to one. These integrals are called "masters" . A priori they have to 
be calculated by independent means. The attempt to eliminate a master usually trivializes the system of 
equations. 

The other integrals in Figure 13 and the seven-propagator configurations found by cancelling a line 
from M 35 , M 36 either have triangles or they are trivial like the last picture in Figure 13. However, 
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according to equation (85) a one-loop subintcgral leaves behind a propagator with non-integer dimension; 
in Figure 14 we marked this by a cross. If this affects one of the a lines of a triangle subgraph, equation 
(84) ceases to be helpful. As a consequence, in the last step a variety of Ti cases with modified propagators 
again have to be attacked by the Laporta idea. One can eventually backsubsitute starting from a general 
result for Ti(l + a\ e, 1 + a 2 e, 1 + a 3 e, 1 + 04 e, 1 + a 5 e) due to [15] and explicit results for the masters 
M 35 , M 36 [16]. 

The integral in the third picture in Figure 12 can be dealt with in the same manner: The bubble 
integral leads to the three-loop NO topology with a cross on one of the outer lines. Once again, the 
Laporta idea is needed to further reduce this. Foreseeably, there is a master integral for which we choose 
NO(e, 1, 1, 1, 1, 1, 1, 1). In Figure 14 we list the more trivial seven-propagator configurations, of which 
only the first case needs Laporta reduction. 





Figure 14 

Finally, topologies B, C, D arc somewhat nicer in that one can always get to the eight-propagator level 
by two successive applications of the triangle rule. They contain the same master integrals M35 , M36 and 
only a BU topology with a cross on the P4 or p 5 propagator (c.f. Figure 3) comes in additionally. 

The recursive procedure is quite laborious because a separate routine is needed for the evaluation 
of every triangle solvable derived topology with a numerator. Further, we had to create tables of M35 
and M36 with two dots and three dot products in the numerator, likewise for the other non-elementary 
cases. These M 35 , M 36 integrals present a formidable computing challenge because the matrix size is of 
the order 1650 * 4400 in polynomials in e. It seems more appealing to attempt a complete numerator 
reduction directly on the topologies A, B, C, D because there are no dots and one needs to consider 
only four distinct mixed products in the numerator so that the matrix results considerably smaller. 

The advantage of our architecture is that the Gaussian elimination need not be pushed to the end; 
it was usually sufficient to eliminate less than half of the integrals in a basis to get definite results for 
any given numerator term. We realised the Gauss algorithm between ANSI C and Mathematica. The 
whole problem addressed in this article could be solved on desktop machines with moderate memory, 
although the total runtime amounted to several weeks. A sample calculation with the programmes will 
be presented in a separate publication, probably on the example of M 36 . 

We conclude this section with a nice observation on the Laporta algorithm: The integrals A, C 
without any numerator are of dimension 1 /q 4 (the fractional part is not important here) and thus power- 
counting finite, and they remain so with numerator pi.q. In our programmes we had not substituted 
explicit results for the master integrals till the very end, but rather kept the rational coefficient function 
that the Gaussian elimination runs up in the corresponding columns. Finiteness of A, C with any of 
the aforementioned numerators constrains the leading orders of the master integrals. The B, D graphs 
without a numerator are divergent due to their triangle subgraphs like e.g. the BE constellation at three 
loops. With numerator p\ or p\ D becomes finite, on B we may use any Pi.q, i G {1 ... 4}. 

The system of equations is strong enough to determine 
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(89) 



in exact agreement with the literature [16] upon conversion to the "G-scheme" . The given orders of the 
e expansion of the masters are also all that our project required. The Laporta algorithm hence turned 
out to be self-sufficient. 
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Finiteness constraints on the level of the nine-propagator integral in the middle of Figure 12 determine 
one order less in M35, M^q 1 and by finiteness of M35, M36 constellations (we may introduce Pi-q/pf on 
any line) we can only fix the first two terms in M 35 while no condition on M 36 ensues. 

7 Conclusions 

In this work we discussed the Konishi anomaly as an operator mixing problem in M = 4 SYM with gauge 
group SU(N). Up to normalisation the two operators that mix are the unintegrated chiral superpotential 
B and Yang-Mills action F of the N — 1 supcrhcld formulation of the theory. Both are not finite on 
their own, but the linear combination O = F — 4 g £> is protected. The second mixture is fixed by the 
conformal properties of two-point functions, i.e. orthogonality to the protected operator and the form of 
its two-point function. We considered the first two non-trivial orders in perturbation theory using SSDR 
as a regulator and multiplicatively renormalised by Z-factors. In the result 

K = Z B B + ^-Z F F (90) 

both renormalisation factors acquire singular higher loop corrections. Further, Zp also receives non- 
vanishing finite corrections at higher loops. Our work fully confirms the general considerations about 
the singularities of the Z-factors put forward in [9] and extends the leading order perturbative analysis 
presented there. The descendant operator K has the anomalous dimension of the Konishi operator as 
required by supersymmetry. 

In conclusion, in our framework the Konishi anomaly is not one-loop exact in contradiction to the 
comment after equation (2.110) in [8]. As an explanation we remark that higher loop mixing coefficients 
are usually scheme dependent. The approach quoted in [8] supposes a definite prescription for the 
renormalisation of the coupling constant appropriate to general J\f = 1 theories, while our scheme is 
tailor-made for the conformal M = 4 case. A general discussion in rigorous perturbation theory wide 
enough to reconcile the two contrary points of view is given in [17]. 

The technically hardest part of the project was to elaborate the four-loop two-point correlator (BF) g s, 
which we achieved by the Laporta algorithm, i.e. integration by parts paired with Gaussian elimination. 
We hope to separately publish the computer programmes developed to this end. Apart from the interest 
inherent to the renormalisation properties of the Konishi anomaly, the four-loop correlator is a vital piece 
of the calculation of the four-loop anomalous dimension of the operator K advocated in [18] along the 
lines of [14]; agreement with the existing result [19] would at the same time confirm the correctness of 
the method (thus the absence of a second "anomaly" in the supersymmetry variation of K) and further 
vindicate the aptitude of the thermodynamic Bethe ansatz [20] to describe "wrapping corrections" in the 
Af = 4 operator spectrum problem, i.e. the regime in which the dilatation generator in the sense of [21] 
becomes inapplicable because the loop order exceeds the spin chain length. 
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Appendix: Feynman rules and conventions 



The N = 4 SYM action formulated in terms of Af = 1 supcrfields has the form 

SV =4 = J d 4 xd 2 6d 2 6Tr(e 9V $ I e- 9V <P I ) (91) 

+ | J d 4 x L d 2 6 e IJK Tr(<t> I [<j ) J ^ K }) + c.c. 
- ^ / d4x L d 2 Tr{W a W a ) 
_ ±Jd*xd 2 ed 2 eTr(^(-±D 2 v)(-^D 2 \ 

+ J d A xd 2 9d 2 9Tr((b + b)L9 V ((c + c) - cothL* v (c - c) )) . 
The definition of the non-abelian field strength multiplet VF Q is 

^ Q = - l -D a D & {f v D a e-<> v ) . (92) 

We choose Fermi-Feynman gauge by putting a = 1: the quadratic part of the YM action becomes 
+1/2 J" VUV. The action has the cubic and quartic YM self-interaction vertices 
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|^ y ^d^^TrQ^^F)]^^^ [v,(D a V)]) 



For completeness, the first two ghost vertices are 

+ | J d 4 xd 2 6d 2 6Tr((b + b) V,{c + c) ) 



|^ / d 4 zd 2 0d 2 0Tr Q(6 + 6), v] [v, (c - c)] ) 



(93) 
(94) 
(95) 

(96) 
(97) 



W.r.t. the definitions in [3], the entire YM and ghost part of the action changes sign, which is related to 
the spinor convention in the integration measure. In addition, V o — V. The ghost propagators are 



(c(2)&(i)) = (6(2) C (i)) = mm) 



(98) 



In Minkowski space the action (91) appears in a weight factor e under the path integral. Instead 
of putting the i onto each vertex we Wick rotate before determining the Feynman rules and proceed in 
positive Euclidean signature. 

The superfields carry an adjoint representation of the gauge group SU(N), and the (Hermitean) 
generators and the structure constants satisfy the relations 



Tr(? ia T'^) S a ^ jibe jabd 2 J\J fi c d 

The spinor convention is: 

V> a = Va = e a p^, ei2 - 1, e 



ae 01 — 



(99) 



(100) 



and exactly the same with dotted indices. Complex conjugation replaces an undottcd by a dotted index 
and vice versa; however, it does not exchange lower and upper position. 

The 2x2 sigma matrices are Hermitian. The a matrix is obtained from a a a by raising of both indices 
as defined by the last equation. They satisfy the following relations: 



(101) 
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The partial spinor derivative obeys 

d a e^ = 5i, d & e^ = si. (102) 

Derivatives with upper indices are defined by the raising rule given in the first equation in (100) and the 
same with dotted indices. (In case of doubt w.r.t. to signs the best strategy is always to put the indices 
into standard position: up on spinors and down on derivatives.) The superspace covariant derivatives are 

D a = d a + i9 a d aa , D a = ~d a - i9 a d aa , d aa = d^G^aa ■ (103) 

We define the square of a spinor without any weight factor: 

e 2 = 9 a e a , e 2 = ej h , (104) 

and similarly for the product of two different spinors. Consequently, 

D.D6 2 = --D 2 9 2 = 1, D.D9 2 = --D 2 9 2 = 1. (105) 
4 4 

Under the x-integral we may thus identify 

d 2 9 = D.D, d 2 9 = D.D . (106) 

SSDR means to modify the dimension of space, but to leave the two-component spinor algebra un- 
touched. In D = 4 — 2 e the basic propagator PL-, (37) becomes 

% = - 4 ^Jh mmh)- (107) 

The propagators carry Kronecker delta symbols for colour and flavour indices, of course. We do not 
usually write these in order to unclutter the notation. 
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